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Abstract 

Where biological datasets are spatially limited, abiotic surrogates have been advocated to inform objective planning for 
Marine Protected Areas. However, this approach assumes close correlation between abiotic and biotic patterns. The Solitary 
Islands IVlarine Park, northern NSW, Australia, currently uses a habitat classification system (HCS) to assist with planning, but 
this is based only on data for reefs. We used Baited Remote Underwater Videos (BRUVs) to survey fish assemblages of 
unconsolidated substrata at different depths, distances from shore, and across an along-shore spatial scale of 10 s of km (2 
transects) to examine how well the HCS works for this dominant habitat. We used multivariate regression modelling to 
examine the importance of these, and other environmental factors (backscatter intensity, fine-scale bathymetric variation 
and rugosity), in structuring fish assemblages. There were significant differences in fish assemblages across depths, distance 
from shore, and over the medium spatial scale of the study: together, these factors generated the optimum model in 
multivariate regression. However, marginal tests suggested that backscatter intensity, which itself is a surrogate for 
sediment type and hardness, might also influence fish assemblages and needs further investigation. Species richness was 
significantly different across all factors: however, total MaxN only differed significantly between locations. This study 
demonstrates that the pre-existing abiotic HCS only partially represents the range of fish assemblages of unconsolidated 
habitats in the region. 
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Introduction 

To adequately represent the entire range of biota present, 
conservation planning for sub-tidal marine ecosystems should be 
based on comprehensive understanding of species, habitats, 
ecosystems and associated ecological processes [1-2]. However, 
in most cases, these data are not available, and hence much of the 
decision-making for planning and management of Marine 
Protected Areas (MPAs) relies on incomplete species inventories 
and distributions [1,3-5]. New approaches that better integrate 
abiotic data into overall planning have consequently been 
advocated [1,4,6-7]. Physical and biophysical surrogates for biotic 
communities are now commonly used in planning of MPAs as they 
are more easily measured, categorised and geo-referenced [8-9]. 
Recent improvements in sonar equipment and techniques now 
provide high-resolution imagery of the physical properties and 
topography of the seabed [10-12]. However, the use of abiotic 
characteristics in planning carries the assumption that biotic 
patterns, which are the end-points of conservation management, 
consistently correlate with abiotic variables. While abiotic 
variables have sometimes been shown to accurately predict biotic 
patterns [13-16], this is not always the case [8]. Testing needs to 



be on a scale relevant to MPA planning [8,16-17], and better 
information on biotic distributions is required for integration into 
existing planning and management of MPAs [4,16]. 

Sub-tidal sediments are the most extensive beiithic habitat 
worldwide [18]. However, there is a paucity of data for fish 
assemblages when compared with adjacent reef habitats [19]. The 
fact that unconsolidated habitats are critical for many commercial 
fisheries worldwide further emphasises the need for greater 
understanding of these habitats and their associated biota [20— 
21]. Often, the use of abiotic surrogates is an unavoidable strategy, 
due to the lack of biotic data [1,3-4,22], and the lack of fish 
assemblage data in unconsoUdated habitats restricts objective 
selection of ecologically important areas for protection. 

Most data on demersal fish assemblages of unconsolidated 
habitats have been collected by trawling [20,23-24], an approach 
that is generally incompatible with MPA objectives [25]. In this 
study, we used Baited Remote Underwater Video (BRUV), which 
is a low-impact, remote method with increasing application to a 
range of sub-tidal habitats [25-33], to assess the effectiveness of a 
Habitat Classification Scheme (HCS) for representing fish 
assemblages of unconsolidated habitats. The current HCS for 
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Figure 1. The location of the sixty-five individual BRUVs deployments in the Solitary Islands Marine Park (SIMP), New South Wales, 
Australia. The northern sampling site is approximately adjacent to Station Creek, and the southern sampling site is approximately adjacent to Bare 
Bluff. The insert box (right side of image) displays the full extent of the SIMP in grey, with the complete sampling area indicated within the red box. 
The second insert box (top of image) shows the study position in relation to Australia. 
doi:10.1371/journal.pone.0096798.g001 



marine parks in New South Wales (NSW), Australia, is an 
important tool for the design and planning of multi-purpose MPAs 
[34] and uses a hierarchical classification, primarily based on 
abiotic factors: habitat type (unconsolidated sediments and hard 
substrata) and depth (shallow: <25 m, intermediate: 25-50 m, 
and deep: >50 m) [16,34]. 

The Solitary Islands Marine Park (SIMP) in the Tweed- 
Moreton Bioregion of northern New Soutii Wales, eastern 
Australia, comprises both State waters (to 3 nautical miles from 
shore), and the adjacent Solitary Islands Marine Reserve (SIMR) 
in Commonwealth waters. The SIMP is unique among NSW 
marine parks in that it includes several islands up to 12 km from 
the mainland coastline. The dominant system of currents, which 
includes strong influence of the East Australian Current (EAC) 
[35-36], cooler coastal counter-currents, and periodic upwelling 
[37], results in an overlap between temperate, tropical and 
subtropical endemic biota in the region [38—39]. The SIMP has 
extensive areas of both rocky reef and unconsolidated sediments 
from intertidal areas to depths >70 m [40]. Comprehensive 
surveys offish associated with hard substrata have detected strong 
depth-related patterns in assemblage structure [16]: depth 
categories for the HCS for the SIMP (shallow: <25 m, interme- 
diate: 25-50 m, deep: >50 m) were consequently refined to reflect 
these patterns. Distance from the mainland coast, which is 
independent of depth due to the presence of shallow fringing reef 
around the islands, also influences patterns of reef fish assemblag- 
es, and HCS therefore incorporates categories to represent these 
patterns (Inshore <1.5 km. Mid-shelf = 1.5—6 km, Offshore > 
6 km). 

Unconsolidated sediments dominate sub-tidal habitats in the 
SIMP [40], but there are currently no data to test whether the 
HCS adequately represents the diversity of fish communities of 
these habitats. However, preliminary BRUVs data, both from an 
examination of the effects of proximity to reef on fish assemblages 
[30], and from haphazard deployment across the SIMP (Malcolm, 
unpublished data), suggest a substantially different suite of species 
to those found on rocky reefs. 

The focus of this research, therefore, was to examine fish 
assemblage structure of unconsolidated habitats to determine 
patterns over factors known to affect their reef counterparts (depth. 



distance from shore), and to determine how consistent they are 
over medium spatial scales (km). As data on the broad composition 
of unconsolidated habitats has recently been generated using 
swath acoustic multi-beam sidescan sonar [40], we also explored 
the relationship between fish assemblages and physical habitat 
characteristics. The specific aims were to: 1) quantify patterns of 
fish assemblages across different depths, distances from shore and 
over a medium spatial scale; and 2) determine which of these, and 
other environmental factors (backscatter intensity, fine-scale 
bathymetric range, and rugosity), are important in structuring 
fish assemblages of unconsolidated habitats in the SIMP using a 
multivariate regression model. 

Methods 

The SIMP covers an area of approximately 71,000 hectares, 
within a broad ecotone between tropical and temperate assem- 
blages [16,40]. The marine park is zoned for multiple use, with 
sections either fully protected (Sanctuary Zone), protected from 
commercial trawling (Habitat Protection Zone), or having no 
protection (General Use). Sixty-five Baited Remote Underwater 
Video (BRUV) deployments were completed within the SIMP 
during the Austral winter of 2011 (Fig. 1, Table 1). The New 
South Wales Marine Park Authority provided the permit for this 
research. As all deployments were within the SIMP boundaries, no 
specific permit was required for this research regarding private 
land access. Similarly, as all sampling was non-extractive and no 
collection of endangered or threatened species or other vertebrate 
fauna was performed, no specific permission was required in this 
context. BRUVs deployments were spread between two broad 
cross-shelf transects ('locations') located in the north and south of 
the SIMP, where the SIMP is longitudinally at its widest, and on 
seafloor that had previously been swath mapped using a 125 kHz 
GeoSwath interferometric sidescan sonar [40]. BRUV deploy- 
ments were also spread across depth and distance from shore 
gradients that provided replication across the existing HCS 
categories. The design was not fully balanced due to constraints 
around areas that had been swath mapped: however, this did not 
preclude statistical modelling. Prior to the field work, coordinates 
for each deployment were uploaded to a handheld GPS to 



Table 1. Factor groups, 


levels and number of BRUVs replicates sampled. 






Factor 


Level 


Number of Replicates 


Depth 


10-30 m 


18 




30-50 m 


24 




5CH-m 


23 


Location 


North 


29 




South 


35 


Backscatter Intensity 


0-50 


50 




50-100 


6 




100-150 


8 


150+ 1 



doi:l 0.1 371 /journal.pone.0096798.t001 
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Table 2. Marginal and sequential test results for the distance based linear modelling (DistLM) procedure using forward selection 
and 4999 permutations, examining four factors of interest within the Adjusted model, and two within the AlC model. 





MARGINAL TESTS 


Selection criterion 


Variable 


SS(trace) 


Pseudo-F 


P 


Prop. 




Adjusted 


Depth 


8086.4 


8.8324 


0.0002 


0.12296 






Location 


4927.5 


5.1026 


0.0006 


7.49E-02 






Backscatter Range 10 m 


3364.4 


3.3966 


0.0092 


5.12E-02 






Rugosity 


1205.7 


1.1766 


0.3168 


1.83E-02 




AlC 


Depth 


8086.4 


8.8324 


0.001 


0.12296 






Location 


4927.5 


5.1026 


0.001 


7.49E-02 






Backscatter Range 10 m 


3364.4 


3.3966 


0.011 


5.12E-02 






Rugosity 


1 205.7 


1.1766 


0.314 


1.83E-02 




SEQUENTIAL TESTS 


Adjusted R^ 


Variable 


SS(trace) 


Pseudo-F 


P 


Prop. (%) 


Cumul. (%) 


0.10904 


Depth 


8086.4 


8.8324 


0.0002 


0.12296 


0.12296 


0.17173 


Location 


4909.8 


5.7687 


0.0002 


7.47E-02 


0.19761 


0.17258 


Rugosity 


904.3 


1 .0636 


0.3836 


1.38E-02 


0.21137 


0.17937 


Backscatter Range 10 m 


1 269.2 


1.5051 


0.1998 


1.93E-02 


0.23066 


AlC 


Variable 


SS(trace) 


Pseudo-F 


P 


Prop. 


Cumul. 


445.24 


Depth 


8086.4 


8.8324 


0.001 


1.23E-01 


0.12296 


441.45 


Location 


4909.8 


5.7687 


0.001 


7.47E-C2 


0.19761 



Significant values are shown in bold. 
doi:l 0.1 371 /journal.pone.0096798.t002 



accurately position BRUVs in the field. There was a minunum 
distance of 200 m between deployments [16,30,41] and all 
deployments were made at least 400 m from the nearest reef, a 
distance that has been shown to be beyond the foraging range of 
most reef-associated fish [30]. 

Field methods 

Each BRUV drop consisted of a mini-DV video camera with a 
wide angle lens, within an underwater housing with flat acrylic end 
ports, an attachment frame, a bait pole and mesh bag with bait, 
and a rope and float system linking the unit to the surface [30,41]. 
Approximately 1 kg of pilchard [Sardimps neopikhardus) was mashed 
into the bait bag, approximately 1.5 m from the lens of the 
camera, to attract fish to the viewing area in front of the camera. 
Each housing and bait pole was bolted to the attachment frame so 
that fish could be viewed in a horizontal orientation to the 
substratum. The field-of-view was standardised to approximately 
2 m behind the bait, to minimise the effects of water clarity on 
measures of relative abundance [30,41]. All files were analysed by 
the same person (AS) to minimise variability in estimating the 
field-of-view. Tapes were converted to digital format (avi) using 
Adobe Premiere Elements (Version 10, Adobe Systems Pty Ltd), at 
a resolution suitable for fish identification (720x576) [30]. 

Video interrogation 

FUes were analysed using Eventmeasure (SeaGIS Pty Ltd, 
Version 3.31), and the identity and MaxN of each fish species was 
recorded. MaxN is the maximum number of fish of each species 
within the field-of-view at any one time during the 30-min 
recording, which removes the possibility of recounting the same 
fish. As counts reflect relative abundance and not density, we 
expect our data to be robust to variability in observing the field-of- 
view [41]. 



Habitat factors 

The placement for all BRUV deployments was determined 
using GIS high-24 resolution imaging software (ArcGIS v9.3, 
ESRI software, USA) which incorporated metrics for the factors of 
primary interest: depth; distance from shore; and along-shore 
location (as distance from latitude -29 S - hereafter termed 
location). Additional factors, including fine-scale bathymetric 
range, rugosity and backscatter were also considered. Bathymetry 
and backscatter data (5 m ceU size) were provided by the NSW 
Office of Environment and Heritage (2012). Rugosity was derived 
from the bathymetry layer using the Benthic Terrain ModeUer 
extension 1.2 [42]. These datasets were interrogated at a 10-m 
radius surrounding the BRUV point location using the Buffer tool 
in ArcGIS (v9.3) and Zonal Statistics tool (ArcGIS Spatial 
Analyst). The Zonal Statistics tool used the 10-m buffer radius 
to calculate ceU statistics (mean, minimum, maximum, range, 
standard deviation and sum) for each layer within this zone. 
Results were exported and incorporated into the variable data for 
subsequent statistical analyses. 

Statistical methods 

Both multivariate and univariate analyses were performed using 
procedures in the PRIMER 6.0 software [43], with the 
PERMANOVA+ add-on [44-45]. Distance-based linear model- 
ling (DistLM) was conducted to model the percentage of overaU 
variation in fish assemblage structure accounted for by each 
environmental factor of interest [44—47]. This procedure utilised a 
permutational approach, whereby the similarity between samples 
was calculated using the Bray-Curtis similarity measure, and a 
forward-stepping selection procedure was applied to regress the 
suite of environmental factors on the resultant matrix. Actual 
values for depth, distance from shore, location, backscatter and 
bathymetric range (within 10-m radius), and rugosity (averaged 
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Figure 2. Distance-based redundancy analysis biplots representing raw Pearson correlations for habitat variables. Vectors are 
overlaid to represent the different environmental variables most important in each modelling approach. Length and direction of vectors indicate the 
strength and direction of the relationship. Numbers 1-4 within the key indicate backscatter intensity. 
doi:1 0.1 371 /journal.pone.0096798.g002 



within 10-m radius) around each BRUV position, provided 
continuous abiotic data for the regression. The factor with the 
highest percentage influence on fish assemblage structure is 
selected first in this procedure, followed by the next highest 
ranked, in sequential order, until no further improvement to the 
model is made. Both adjusted R" and Akaike Information 
Criterion (AIC) [48] selection criterion were tested in the analysis 
to provide a comprehensive evaluation of appropriate predictor 



factors to include in the model. As regression-based models are 
sensitive to correlation between factors [49] , environmental factors 
with a correlation coefficient of ^0.7 were noted, and one of the 
correlated pair was excluded from the analysis [49-50]. Such 
correlations were found between depth and distance from shore, 
and between bathymetric range and rugosity, so distance from 
shore and bathymetric range were discarded from the analysis, 
and depth and rugosity, along with location and backscatter, were 
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Figure 3. Distance-based redundancy analysis bipiots representing raw Pearson correlations for fish species. This is as seen in Fig. 2, 
however here vectors for the 13 most influential fish species in the analysis are overlaid. Length and direction of vectors indicate the strength and 
direction of the relationship. Numbers 1-4 within the key indicate backscatter intensity. 
doi:10.1371/journal.pone.0096798.g003 



retained. These variables were retained as they were considered to 
be the most intuitive and interpretable of the measures available. 

Distance based redundancy analysis (dbRDA) bipiots were 
generated to visually display the direction and magnitude of the 
relationship between habitat factors and individual fish species 
[47]. If DistLM indicated a significant depth effect, BRUV 
positions were categorised for that factor using the existing SIMP 
habitat classification system (shallow = 10-30 m, intermediate 



= 30-50, deep = 50+). If location was significant, BRUV positions 
were categorised as north or south. If backscatter was significant it 
was analysed by categorising intensity range (maximum value 
minus minimum value) into four levels (0-50, 50-100, 100-150, 
and 150+) as this could be considered most analogous to sediment 
complexity from the values available. Likewise, if rugosity was 
significant it was analysed by categorising intensity range into four 
categories (0.00001-0.0001, 0.00011-0.001, 0.00101-0.01 and 
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Table 3. Family, genus, species, common name, and Pearson correlations of environmental variables with each of the four dbRDA 
axes, for both the adjusted R2 and AlC models. 




Family 


Genus 


Species 


Common Name 


Correlation 




RDA1 








Adjusted R^ 


AlC 


Myliobatididae 


Myliobatis 


austral is 


Eagle Ray 


0.44 


0.55 


Platycephalidae 


Platycephalus 


longispinis 


Long-Spine Flathead 


0.31 


0.33 


Scombridae 


Sardo 


austral is 


Australian Bonito 


0.26 


0.22 


Urolophidae 


Trygonoptera 


testacea 


Common Stingaree 


0.19 


0.30 


Rachycentridae 


Rachycentron 


canadum 


Black Kingfish 


0.19 


0.32 


Platycephalidae 


Platycephalus 


caeruleopunctatus 


Blue-Spotted Flathead 


-0.54 


-0.47 


Triakidae 


Mustelus 


antarcticus 


Gummy Shark 


-0.25 


-0.25 


Rhinobatidae 


Aptychotrema 


rostrata 


Eastern Shovelnose Ray 


-0.25 


-0.15 


Aracanidae 


Anoplacapros 


inermis 


Eastern Smooth Boxfish 


-0.04 


-0.11 


Tetraodontidae 


Lagocephalus 


sceleratus 


Silver Pufferfish 


-0.10 


-0.16 


RDA2 


Tetraodontidae 


Lagocephalus 


sceleratus 


Silver Pufferfish 


0.24 


0.26 


Aracanidae 


Anoplacapros 


inermis 


Eastern Smooth Boxfish 


0.21 


0.19 


Pinguipedidae 


Parapercis 


stricticeps 


White-Streaked Grubfish 


0.12 


0.08 


Scorpaenidae 


Centropogon 


austral is 


Fortesque 


0.03 


0.08 


Urolophidae 


Trygonoptera 


testacea 


Common Stingaree 


-0.34 


-0.16 


Rachycentridae 


Rachycentron 


canadum 


Black Kingfish 


-0.34 


-0.10 


Sillaginidae 


Sillago 


spp. 


Whiting spp. 


-0.34 


-0.29 


Myliobatididae 


Myliobatis 


austral is 


Eagle Ray 


-0.31 


-0.10 


Rhinobatidae 


Aptychotrema 


rostrata 


Eastern Shovelnose Ray 


-0.19 


-0.20 


Rhinobatidae 


Trygonorrhina 


fasciata 


Eastern Fiddler Ray 


-0.19 


-0.16 


RDA3 


Rachycentridae 


Rachycentron 


canadum 


Black Kingfish 


0.27 




Urolophidae 


Trygonoptera 


testacea 


Common Stingaree 


0.17 




Tetraodontidae 


Lagocephalus 


sceleratus 


Silver Pufferfish 


0.16 




Aracanidae 


Anoplacapros 


inermis 


Eastern Smooth Boxfish 


0.09 




Paralichthyidae 


Pseudorhombus 


arsius 


Large-Tooth Flounder 


-0.20 




Triakidae 


Mustelus 


antarcticus 


Gummy Shark 


-0.16 




Platycephalidae 


Platycephalus 


caeruleopunctatus 


Blue-Spotted Flathead 


-0.13 




Pinguipedidae 


Parapercis 


stricticeps 


White-Streaked Grubfish 


-0.12 




RDA4 


Paralichthyidae 


Pseudorhombus 


arsius 


Large-Tooth Flounder 


0.39 




Myliobatididae 


Myliobatis 


austral is 


Eagle Ray 


0.17 




Rhinobatidae 


Aptychotrema 


rostrata 


Eastern Shovelnose Ray 


0.10 




Pinguipedidae 


Parapercis 


stricticeps 


White-Streaked Grubfish 


0.03 




Platycephalidae 


Platycephalus 


caeruleopunctatus 


Blue-Spotted Flathead 


-0.21 




Urolophidae 


Trygonoptera 


testacea 


Common Stingaree 


-0.10 




Aracanidae 


Anoplacapros 


inermis 


Eastern Smooth Boxfish 


-0.08 




Triakidae 


Mustelus 


antarcticus 


Gummy Shark 


-0.08 





doi:l 0.1 371/journal.pone.0096798.t003 



0.01001+) as it could be considered most analogous to benthic 
topographic complexity. These biplots were used to visually assess 
the relationship between existing categories and significant factors. 

Univariate permutational multivariate analysis of variance 
(PERMANOVA) (using Euchdean distance as the similarity 
measure) was performed to test for significant differences in 
Species Richness and Total MaxN between factors that were 
significant within the DistLM regression. Species Richness is 



calculated as the total number of different species viewed on each 
replicate video, and Total MaxN is calculated as the sum of MaxN 
for all species per replicate. Differences across bathymetry and 
rugosity categories were not analysed as the range of values was 
very low. 

One-way PERMANOVA [44-45] was used to test for statistical 
differences in assemblage structure for those factors that were 
significant in the DistLM regression. Each test was run using 4999 



PLOS ONE I www.plosone.org 



7 



May 2014 | Volume 9 | Issue 5 | e96798 



Spatial Patterns of Sedimentary Fish Assemblages 



20 



'15 



X 



Bio 



5 

CO 
CO 
CD 

o 

a: 

Q) 

'o 

Q.2 
CO 



Depth 



i 



X 



i 



X 



Backscatter Intensity 



10-30m 



30-50m 



Depth 



50+m 



0-50 50-100 100-150 150+ 

Backscatter Intensity 



Figure 4. Mean species richness and total MaxN (±SE) at three depths and four different levels of backscatter intensity. The graph for 
depth also show differences between northern (black bars) and southern (grey bars) locations within these factors. Data from northern and southern 
locations has been combined for backscatter intensity due to unbalanced numbers of replicates within the four levels of intensity. 
doi:1 0.1 371 /journal.pone.0096798.g004 



permutations. Where significant results were generated, post-hoc 
pair-wise contrasts were performed to establish where differences 
in assemblage structure were occurring. 

Multivariate assemblage structure was visually examined using 
non-metric multi-dimensional scaling (nMDS) ordination. Data 
were square-root transformed, and a dummy variable was added, 
prior to the generation of Bray-Curtis similarities. The resultant 
nMDS was used to examine differences in assemblage structure 
across factors that were significant in the DistLM analysis. Depth- 
related patterns were visually examined to see how well they fitted 
the current depth-based habitat classification categories [16]. 

Results 

General 

A total of 16 teleost and elasmobranch species, from 12 families, 
was recorded from BRUVs deployments. Platycephalidae was the 



most speciose family (3 species). Three species were numerically 
abundant and ubiquitous across all BRUV positions in the study - 
Platycephalus caeruleopunctatus, P. longispinis and Aptychotrema rostmta. 
Schooling Sillago spp. were numerically abundant in some 
positions, and recorded across all depths, but were highly variable 
between positions and between depths. 

Correlations between habitat variables and assemblage 
structure 

Marginal tests in distance-based linear modelling (DistLM) 
showed that, of the four environmental factors retained for 
analysis, three were significant in isolation (Table 2) for both the 
adjusted R and AIC selection criteria. However, the optimum 
model for adjusted R^ included all four factors: for AIC, only 
depth and location were included in the optimum model. The 
adjusted R^ model explained 23.1% of the total variation, while 
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Table 4, Summary of one-way PERMANOVA results for the ana 
differences in assemblage structure across the different factors. 
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the AIG model explained 19.8% (Fig. 2, Table 2). The first two 
dbRDA axes explained 99.6% of the fitted variation for the 
adjusted model, and 100% for the AIG model. In the context of 
hypothesis generation, the more parsimonious AIC model does 
not allow for further exploration of factors outside depth and 
location, thus it was worthwhile to further explore the adjusted 
model despite it being less parsimonious. Backscatter was found to 
be significant in the marginal tests despite it being highly 
unbalanced across the study design, and thus also warranted 
further exploration. Raw Pearson correlations of each factor with 
each dbRDA axis for the adjusted R^ model showed depth (p = 
— 0.90) correlated closely with the first dbRDA axis, with 
backscatter intensity also partially correlated (p = 0.62). The 
second axis was a combination of location (p = —0.85) and depth 
(p=-0.42). For the AIC model depth (p=-0.99) was almost 
perfectly correlated with the first axis, while location (p = —0.99) 
was similarly correlated with the second axis. 

By examining the raw Pearson correlations for different species, 
we could identify those most responsible for driving the 
assemblage response to each environmental factor: these are 
displayed as vectors in Fig. 3. The data for the adjusted R^ model 
must be interpreted somewhat cautiously, as the first dbRDA axis 



represents a combination of depth and backscatter, while the 
second axis is a combination of location and depth. Therefore, 
determination of the factor primarily responsible for the distribu- 
tion of each species is partially confounded. Some of the species 
and their respective correlation values for both the adjusted R^ 
and AIC models are displayed in Table 3. Species with the four 
most positive and four most negative correlations with each axis 
are displayed, as these are the species most responsible for driving 
patterns in the assemblage. 

Univariate measures 

Species richness was generally higher in the south of the park 
(Fig. 4). The highest species richness by depth categories was for 
shallow deployments in the south of the park and the second 
highest was for deep deployments in the south. As the number of 
replicates was unbalanced across backscatter intensity levels 
(Table 1), all deployments were pooled for graphical representa- 
tion. The highest species richness occurred in the 50-100 category 
with the lowest being in the 150+ category. 

PERMANOVA revealed a significant difference between 
species richness for all four factors of interest (Table 4). Pairwise 
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Table 5. Results of post-hoc pair-wise contrasts of species richness and total MaxN, and assemblage structure (PERMANOVA) for 
each pair of levels across the factors. 
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contrasts revealed significant differences between shallow and 
intermediate, and intermediate and deep categories (Table 5). The 
only significant contrasts for backscatter were between levels 0-50 
and 50-100 (Table 5). 

Total MaxN values were slightly higher in the south of the park, 
and showed trends similar to those seen for species richness (Fig. 4). 
The highest Total MaxN was recorded in the shallow southern 
category and the second highest in the deep southern category. As 
with species richness, highest Total MaxN across backscatter 
categories was for 50-100. Of the three factors of interest, 
PERMANOVA found a significant eflFect only for the latitude 
(Table 4). 

Assemblage structure 

PERMANOVA revealed significant effects for all three selected 
factors (Table 4). Pairwise contrasts between different levels of 
depth were all significant. Pairwise contrasts across backscatter 
categories were complicated by unbalanced numbers of samples in 
each level, generating a large number of unique permutations 
between some pairs of samples, and relatively few between other 
pairs (Table 4), making it impossible to make statistical inferences 
at a level £0.05 for several pairs. 

Some trends were evident in the iiMDS ordination (Fig. 5) with 
the majority of shallow inshore and mid-shelf sites grouping in the 
top half of the plot, and the majority of the deep offshore sites in 
the lower half In addition, there was some along-shore separation 



between assemblages, with those in northern locations tending to 
the left of the ordination, and those from southern locations 
tending to the right (Fig. 5). 

Discussion 

Patterns of fish assemblage structure in unconsolidated habitats 
within the SIMP are clearly influenced both by depth and factors 
operating over medium spatial scales (i.e. the distances between 
our locations). Additionally, backscatter intensity, while not a 
primary driver, was highlighted as potentially explaining some of 
the diflerences in fish assemblages across the scales of the study. 

Assemblage patterns and correlations with 
environmental variables 

Worldwide, fish assemblage patterns of both reef and uncon- 
solidated habitats have often been found to be correlated with 
depth [16,21,23,28,32,51-52] and distance from shore [26,34,53], 
and these factors have consequeiitiy been used as biophysical 
surrogates for conservation planning in MPAs [16,34]. Fishes are 
an important ecological and socio-economic component of marine 
diversity in their own right [54] but can also reflect the distribution 
of benthic communities [16,34]. Distinct patterns in the structure 
of fish assemblages of unconsolidated habitats have also been 
documented in previous research in coastal and shelf waters of 
eastern Australia [14,19-20]. For example, [19] found distinct 
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Figure 5. Non-metric multi-dimensional scaling (nMDS) ordination showing the relationship among fish assemblages from three 
depths, three distances from shore, and two locations. Data were square-root transformed prior to analysis. Lines represent 60% similarity. 
Numbers 1-4 within the key indicate backscatter intensity. 
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assemblages in three depth categories (20-30 m, 40-50 m and GO- 
TO m) on unconsolidated sediments offshore from Sydney. They 
also showed that: while broad patterns were evident, none of the 
abundant species consistently showed a preference for one depth 
category; species that were unique to a single depth category were 
found in low abundance (<3); and the numerically dominant 
species were found at all depths [19]. Similarly, [14] concluded 
that assemblage patterns were being driven mostiy by differences 
in the abundance of common species that occurred across multiple 
habitats and depths. Our data show a similar pattern in that the 
species most responsible for differences in assemblage structure are 
also numerically abundant, and were recorded at all depths and 
within different locations in the SIMP. 

The importance of depth in structuring assemblages reflects 
what is known from other habitats in the SIMP and elsewhere. 
Strong depth-related patterns have been detected in reef fishes 
[34], for motile invertebrates of natural [55] and artificial [56] 
substrata, and for habitat-forming benthos [57]. The shallow/ 
intermediate boundary in the current HCS in the SIMP was 
defined at 25 m due to the disjunct patterns of reef fish and 
general benthic assemblages (from coral to sponge-dominated) 
below this depth [58]. 

With strong co-linearity between depth and distance from shore 
within our study, it is difficult to distinguish which factor primarily 
drives the observed patterns. Shelf position is included in the 
current HCS for the SIMP, as distinct cross-shelf patterns have 
been demonstrated for corals [59], molluscs [60] and reef fish [34]. 
Due to the broad morphology of the continental shelf region in the 
SIMP, and the mobile nature of unconsolidated habitats, this 
correlation between factors is not unexpected, and from a practical 
perspective, depth could be considered a proxy for distance from 
shore, or vice-versa, for management purposes. 

The correlation with location in our modelling contradicts 
findings for reef fish communities in the SIMP [16,34] which 
showed no noticeable pattern between assemblage structure over 
similar spatial scales, and it is currently unclear as to the reasons 



for this variation. While broad-scale latitudinal effects have been 
observed in fish communities of unconsolidated sediments 
[24,26,61-62], these patterns are generally evident over a much 
larger spatial scale than in our study (e.g. [26]). It is likely that this 
variation is due to other factors not identified within this study, 
rather than actual latitudinal effects. 

Our modelling also suggests the possibility that sediment grain 
size, for which backscatter intensity is an appropriate surrogate, 
may influence fish communities of unconsolidated habitats. Broad- 
scale bathymetric mapping of the seafloor within the SIMP has 
shown that unconsolidated habitats show distinct variability in 
sediment type, as revealed in acoustic backscatter [40] . Some areas 
show a darker, 'reef-like' backscatter, and are likely to contain 
varying amounts of pebbles and cobbles [40] . These patches of 
denser backscatter, while extensive in some sections of the SIMP, 
were not the primary focus of our study. During this study, 4 
species {Pseudorhombus arsius, Parapercis stricticeps, Centropogon australis, 
Ratabulus diversidens) were found only at sites with a backscatter 
index S50, despite a low number of replicates with tiiis index or 
greater (15) relative to the total (65, see Table 1). Additionally, 
while our sampling design did not consider backscatter intensity as 
an a priori factor, the marginal tests from the DistLM analysis 
generated a sigiiiflcant value for this factor. In the context of 
hypothesis generation, we recommended that habitat types should 
be examined at a finer scale to comprehensively assess the role of 
backscatter intensity (as a proxy for sediment type) in the 
distribution of different species. With obvious and extensive 
differences in sediment characteristics, it is likely that, in addition 
to the other factors identified in this study, sediment type and 
habitat complexity are important factors structuring fish assem- 
blages. Recent research in the SIMP [30] has demonstrated that 
the presence of reef adjacent to sedimentary habitats, even at the 
scale of hundreds of metres, can have a substantial effect on fish 
assemblage structure. It is possible that gravels may act as a 'reef 
surrogate' for some species. 
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When compared to rcifs, unconsolidated habitats, particularly 
in a high energy environment, have lower topographical 
complexity, support fewer sessile biota, and provide limited 
foraging opportunities and protection from predation [63-65]. 
The low species richness observed in this study reflects this, with 
the species complement comprising ~3% of the overall teleost and 
elasmobranch diversity recorded from the SIMP [34] . Despite this 
low species richness, statistical differences were evident over levels 
of each of the four factors we tested, and rare species were mosdy 
restricted to single levels within each factor. In contrast, the only 
significant factor for total MaxN was location, with higher overall 
MaxN values at the southern location. 

The current Habitat Classification System (HCS) for the SIMP 
has been developed and adjusted primarily based on fish 
assemblage patterns associated with hard substrata [16,34]. This 
HCS has been used to systematically examine representation of 
habitat categories as a surrogate for biotic pattern [66-68], under 
existing and different marine park zoning scenarios [69]. In the 
absence of specific data, the same categories have been used for 
conservation planning for sediment-associated biota. There is 
value, both ecologically and from a plarming perspective, in 
recognising that broad assemblage patterns in fish assemblages of 
unconsolidated substrata do exist in the SIMP and require 
consideration for planning elsewhere. That these, at least in part, 
align with the current habitat categories used within this MPA is 
opportune. However, it also indicates that some further refinement 
of the current HCS in the SIMP is required to improve future 
systematic planning. Use of quantitative and systematic planning 
in an MPA, based on suitable and available datasets, is more-likely 
to achieve MPA objectives, as well as address some of the essential 
criteria necessary to effectively evaluate and select areas for 
protection [70-71]. It is also more likely to provide cost-effective 
solutions in terms of area required and/or social impacts [71-72]. 
Biotic data are often spatially constrained, whereas physical/ 
environmental data are often less so [66] and, therefore, 
integrating these can achieve better planning outcomes [73]. 
Without testing at the appropriate scale, however, it cannot be 
assumed that the patterns between fish and physical factors 
observed in this study will occur in other MPAs. Our study 
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